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INTRODUCTION 


The  investigation  of  digital  signal  processing  algorithms  requires 
the  computer  simulation  of  sensor  signals  having  certain  stochastic 
properties.  The  construction  of  such  signals  by  Monte-Carlo  methods 
generally  entails  the  use  of  computer  generated  sequences  of  random 
numbers.  Although  most  computer  library  systems  have  subroutines 
available  for  generating  random  sequences  according  to  a  variety  of 
probability  amplitude  distributions,  almost  without  exception,  their 
elements  are  uncorrelated.  This  enables  the  rapid  and  economical 
simulation  of  combinations  of  primary  and  intefering  signals  (e.g., 
noise)  in  which  the  former  is  deterministic  and  the  latter,  random  and 
uncorrelated.  Alternatively  both  the  primary  and  interfering  signals 
can  be  represented  by  uncorrelated  random  sequences.  An  example  would 
be  a  radar  echo  in  which  the  signal  of  interest  is  random  due  to  random 
fluctuations  in  the  target  cross  section  and  the  interfering  signal  is 
white  Gaussian  noise. 

Between  these  two  extremes  is  a  class  of  interesting  signals  which 
must  be  represented  by  random,  but  correlated  sequences.  This  class 
poses  a  problem  for  computer  signal  simulations  because  in  general  there 
is  no  computationally  feasible  way  to  construct  a  sequence  of  random 
numbers  having  both  a  specified  probability  amplitude  distribution  and  a 
specified  autocorrelation  function.  It  is  not  difficult  however  to 


generate  a  normally  distributed  sequence  with  a  specified  auto¬ 


correlation  function 


BACKGROUND 


The  present  investigation  arose  in  connection  with  an  attempt  to 
simulate  radar  clutter  echoes.  Such  echoes  are  random  and  usually 
correlated  over  several  radar  pulses.  Among  the  detector  envelope 
probability  amplitude  distributions  which  have  been  used  in  radar 
clutter  modeling  and  analysis  are  the  log-normal  and  Weibull 
distributions  [1,2} . 

That  clutter  echoes  can  often  be  described  by  a  log-normal 
distribution  leads  to  a  convenient  way  of  simulating  them.  Because  of 
the  very  simple  transformation  which  relates  a  log-normally  distributed 
random  variable  to  one  which  is  merely  normally  distributed,  it  is  easy 
to  construct  log-normal  sequences  with  specified  autocorrelation 
f unct ions . 

This  problem  has  in  fact  been  addressed  in  the  literature  [3,  4]. 

In  both  of  those  references,  the  method  involved  the  direct  generation 
of  a  normally  distributed  sequence  with  a  specified  covariance  matrix 
followed  by  an  appropriate  transformation.  This  paper  adopts  a  slightly 
different  approach  based  on  a  Discrete  Fourier  Transform  of  the  appro¬ 
priate  spectral  density.  Since  this  can  be  implemented  by  means  of  a 
standard  FFT  routine  it  may  be  faster  than  the  direct  method. 


The  remainder  of  this  paper  will  not  deal  in  detail  with  properties 
of  radar  clutter  per  se,  hut  rather  with  the  techniques  of  generating 
correlated  log-normal  sequences  by  using  correlated  normal  sequences. 
Further  discussion  with  particular  quantitative  application  to  radar 
clutter  modeling  will  be  presented  in  a  forthcoming  paper. 

Tn  the  following  sections  fhe  relationships  between  the  means, 
variances  and  autocorrelation  functions  of  normal  and  log-normal  random 
variables  are  given.  The  generation  of  correlated  normal  sequences 
using  the  power  spectral  density  is  discussed.  The  relationship  between 
the  autocorrelation  functions  is  then  invoked  to  produce  appropriately 
correlated  log-normal  sequences.  Illustrations  of  the  method  are 
presented  in  the  final  section.  The  discussion  is  limited  to  sequences 
which  represent  wide-sense  stationary  processes.  These  are  processes 
whose  means  are  independent  of  time  and  whose  correlation  functions 
depend  only  on  the  time  interval  between  samples* 


MATHEMATICAL  TECHNIQUES 


LOG-NORMAL/ NORMAL  RELATIONSHIPS 


Let  z^  =  z(t)  and  z(t+r)  be  mean  zero  normally  distributed 

random  variables  from  a  wide  sense  stationary  process  {z(t)}  with 
o 

variance  oz  and  correlation  coefficient  p  -p(t).  The  probability 
density  function  of  z^,Z2  Is 


pz(zlz2) 


exp  - 


2'1-p2),2  ['* 


2Pzlz2  + 


(1) 


The  autocorrelation  function  is 

z~zT  =  z(t)z(t+T)  =  R  <t)  -  o2  p (t ) 
*■  z  z  z 


(2) 


in  which  the  over-bar  denotes  statistical  expectation. 


Upon  setting  Iny^  =  z^+a  and  iny2  “  Z2+a  (a  «  constant)  the 
variables  y^  =  y(t)  and  y2  =  y(t+r)  are  found  to  be  log-normally 
distributed  with  joint  density  function: 


py<y1.y2) 


1 

2  /  2 

2i»oz  /l-p  y]y2 

l  f  2 

e*P  2 — 2  C^ny  -e)  -2p(*ny  -a)(iny  -a) 

2(l-pZ)a‘  >-1  12 


+  (iny2-a) 


2] 


(3) 
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By  performing  suitable  integrations  it  is  easy  to  show  that  y^  and  y2 
are  each  separately  log-normal.  Furthermore,  the  statistics  of  the 
{ y}  process  are  simply  related  to  those  of  the  { z)  process  as 
follows: 


(*wL) 


yx  =  y2  =  y  =■  e 


(a) 


_2  _2  _2  2(a4o2) 

yt  -  y2  =■  y  -  e 


2  -2  .-.2  2.pz  °\1 

°y  -  y  -  (y)  -  e  l_e  -  e  J 


(b) 


(c) 


(4) 


2  —  2  —  2  Rz^T  ^ 

Ry(0  -  a/y<T>  +  (y)  -  (y)  e 


(d) 


or 


inR  (x) 

r 


"  ,-.2 

(y) 


in 


(y)  y 


(e) 


Equations  (4,  a-e)  contain  the  key  to  the  procedure.  The  relevant 
statistics  are  related  by  simple  exponential  or  logarithmic 
transformations.  For  example,  if  {y}  is  a  log-normal  process  with 
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autocorrelation  function  Ry(r),  then  under  the  transformat  ion 
any  ■+  z+a  ,  {z}  will  be  a  normal  process  with  autocorrelation 

function  R2(t)  given  by  (4e). 

DISCRETE  FOURIER  TRANSFORMS 

The  primary  tool  required  for  the  generation  of  correlated 
sequences  is  the  finite  Discrete  Fourier  Transform  (DFT).  The 
properties  of  this  transform  follow  logically  from  those  of  the  infinite 
interval  continuous  Fourier  Transform.  The  latter  are  well  known.  A 
few  definitions  are  given  here  mainly  to  establish  the  notation. 

Assuming  the  existence  of  all  integrals,  the  standard  Fourier 
transform  pair  {x(t),  X(ui)}  is 

oo 

*(t)  »  j  X(u))eiu,tdu  (a) 


(5) 


X(u>) 


/  x(t)e  iu)tdt 

—CO 


(b) 


One  important  property  of  the  pair  (5a,  5b)  is  that  for  a  real  time 
function  x(t) 


X(-w)  =  X  (o)>  (real  x(t)) 


(6) 


* 
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where  the  asterisk  denotes  complex  conjugation.  If  in  addition,  x(t) 
is  even  then  X(w)  is  also  real  and  even: 


X(-w)  =  X(w)  (real,  even  x(t)) 


(7) 


Given  a  wide  sense  stationary  process  with  autocorrelation  function 
R(t)  ,  the  spectral  density  S(w)  is  obtained  from  R(t)  by  a  transform 
like  (5b): 

oo 

S(u>)  -  /  (R(T)e"lu)rdT  (8) 

— OO 

and,  since  R(t)  is  real  and  even  so  is  S(m)  . 

One  way  of  effecting  passage  to  the  finite  DFT  is  to  approximate 
(5b)  by  a  discrete  sura  over  an  interval  (-T/2,  T/2)  .  One  sets 
T  “  NAt  and  t  »  nAt  ,  n  “-N/2  . N/2-1  and  evaluates  the 

resulting  summation  at  the  discrete  set  of  frequencies 

a),  =*  kAw  =■  •  The  result  is: 

k  T 


N 

2~1  -12tt  kn 

X(kAui)  ■  At  x(nAt)e  N  (9) 

n*T 
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Standard  practice  Is  to  shift  the  Index  n  so  that  It  runs  from  0  to 
N-l.  Setting  X(kAu>)/At  =  X^  and  x(nAt)  “  xn  in  (9)  and  performing 
this  shift  leads  to  the  usual  form  of  the  DFT  and  its  inverse: 


*k  - 


N-l 

X 

u=o 


-2nikn 


x  e 
n 


N 


(a) 


x 

n 


1_ 

N 


N-l 

l 

k=o 


*k 


2irikn 


(10) 


(b) 


Corresponding  to  (6)  and  (7)  for  real,  and  real  and  even  x(t)  are 
analagous  relations  for 

Xk  *  XN-k  k  =  1,2  f  "  1  (real  xn>  (H) 

Xk  3  XN-k  k  =  i*2'*’*  f  “  1  (real, even  x^)  (12) 

Equations  (10a)  and  (10b)  together  with  the  condition  (12)  are  used 
extensively  in  the  following  paragraphs. 
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GENERATING  CORRELATED  SEQUENCES 


To  generate  a  correlated  sequence  with  a  specified  autocorrelation 
function,  the  first  step  Is  to  compute  the  corresponding  spectral 
density.  For  this  purpose  equation  (8)  can  be  used  provided  the 
Integration  can  be  performed.  Otherwise,  an  approximation  using  (10a) 
with  xn  “  Rj!  =  R(nAx)  will  suffice.  In  either  case,  samples 
Sk  »  S(kAu>)  of  the  spectral  density  will  be  related  to  the  given 
autocorrelation  function  by 


1  N 
k=*o 


2irlkn 

N 


(13) 


which  Is  simply  equation  (10b)  applied  to  the  transform  pair 
With  these  spectral  density  samples  in  hand  generate  two  Independent 
sequences  of  zero  mean  uncorrelated  random  variables  {X^}  ,  {Y^} 
according  to  any  arbitrary  probability  amplitude  distribution,  and  form 
the  complex  random  variable  as  follows: 


Uk  *  Xk  k~I 

Uk  -  Xk  +  1  Yk  k  -  1,2. . .N/2-1 


(U) 


The  Xy  and  Yk  must  be  scaled  so  that  for  each  k 
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(15) 


Using  samples  of  the  spectral  density  obtained  from  the  given 

autocorrelation  function  Rj^  form  a  new  set  of  complex  variables  Zk 
by  scaling  the  this  way: 


ZR  =  k-0,1. .  .N/2 


(16) 


For  k  =  N/2+1,  N/2+2 , • • •  N-l  set 


Z 


k 


(17) 


Now  perform  an  inverse  DFT  on  the  Zk  normalized  as  shown: 


z 

n 


✓  N 


N-l 

1  „  _  N 

N  . Zke 
k«0 


(18) 


The  sequence  (zn)  resulting  from  this  procedure  is  real,  and 
furthermore  has  the  prescribed  autocorrelation  funtion  Rn  .  That  the 
zn  are  real,  follows  immediately  from  (11)  in  light  of  condition  (17). 


-11- 


That  Rn  (corresponding  to  S^)  is  indeed  the  autocorrelation  function 
of  the  zn  can  also  be  easily  demonstrated.  For  the  correlation 
function  of  this  sequence  at  a  lag  i  say  is 


2nkni  -  2tt  im( n-Hl ) 

-  1  _  -=-=*  N  N 

zz,*—  2  2LZe  e 

n  n+l  N  .  it  m 
k,m 

2v  in(  k-m)  -2rc  Imi 

-i  I  V*  eN  e  N 

N  .  k  m 
k  ,ni 


But  Z,  Z  ■  / S,  S  U,  U4  "  / S,  S  5,  ,  since  the  U,  are  uncorrelated 

It  m  It  m  km  km  km  k 

zero  mean  random  variables.  Therefore 


*nznH 


-  E 

“  K.. 


'V.  6 


2tt  lq(k-m)  —  2tt  iml 
N  N 


km 


or 


ZnZn+l  N  \ 
k 


—  2tt  ikJt 
N 


(19) 


In  other  words  z  z  ,  „  is  identical  with  the  inverse  transform  of  S,  , 
n  n  H  k 

and  this  is  precisely  . 
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The  only  remaining  question  Is  the  probability  amplitude 


distribution  of  (zn)  •  Recall  that  {X^}  In  frequency  space  which 
spawned  (zn)  were  zero  mean  uncorrelated  sequences.  But  their 
probability  distribution  was  not  specified.  The  zR  however  are  linear 
combinations  of  the  X^  and  Y^  and  hence  the  Central  Limit  Theorum 
Implies  that  for  some  value  N  ,  of  the  size  of  the  sequence,  {zn}  will 
tend  toward  a  normal  distribution.  The  mean  of  this  distribution  will 
be  zero  and  Its  variance  will  be  given  by  equation  (19)  evaluated  at 
zero  lag: 


o 


2 

z 


1 

N 


N-l 

k-0 


(20) 


An  example  of  a  sequence  {zR}  with  a  given  correlation  function, 
constructed  from  two  uniformly  distributed  sequences  {X^}  and  {Y^} 
will  be  given  in  the  next  section. 

The  technique  described  above  for  generating  correlated  sequences 
together  with  equation  (4)  relating  the  statistics  of  the  log-normal  and 
normal  distributions  enables  the  generation  of  correlated  log-normal 
sequences  with  little  difficulty.  The  method  is  summarized  below. 

a  Given  a  desired  log-normal  process  {y(t)}  with  mean  y  , 
variance  a*  and  autocorrelation  function  Ry(t),  use 
equation  (4e)  to  find  the  autocorrelation  function 


R2(t)  of  the  underlying  zero-mean  normal  process  (z(t)}. 


•  Either  by  direct  Integration  of  equation  (8)  or  use  of  the 
discrete  transform  (10a)  compute  samples  S2(k)  of  the 
power  spectral  density  of  the  {z(t)}  process. 

•  Apply  equations  (14)  -  (18)  to  obtain  the  sequence  (zn). 

In  particular,  generate  the  sequences  {X^}  and  {Yfc}  as 
zero  mean  NORMAL  processes.  Since  a  linear  combination  of 
normally  distributed  random  variables  is  always  Itself 
normally  distributed,  this  guarantees  that  the  zn  are 
normally  distributed  zero-mean  random  variables. 

Z  "^3 

•  Set  yn  -  e  n  for  n  ■  0,1,...  N-l.  The  constant  a  is 
determined  by  y  and  o^.  According  to  equations  (1)  - 
(4)  the  yn  will  be  log-normally  distributed  with  the 
statistics  y  ,  o*  and  Ry(i)  .  The  relationships 
between  the  z  and  y  statistics  are  as  expressed  in 
(4). 

An  example  of  this  procedure  is  given  in  the  next  section. 
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EXAMPLES  OF  CORRELATED  SEQUENCES 


GENERAL 

Two  examples  of  the  preceding  methods  are  presented  in  this 
section.  The  first  one,  Intended  for  comparison  purposes.  Illustrates 
the  generation  of  a  random  sequence  with  a  specified  autocorrelation 
function  but  with  no  attempt  to  shape  the  final  probability  amplitude 
distribution  function.  The  second  illustrates  the  generation  of  a  log¬ 
normal  sequence  with  this  same  autocorrelation.  To  make  the  comparison 
complete,  both  the  mean  and  variance  of  each  of  the  final  distributions 
were  set  equal  to  l. 

To  obtain  reasonable  statistical  samples  an  ensemble  of  100 
sequences  was  generated  for  each  example.  Each  of  the  sequences  had  64 
elements  leading  to  a  total  of  6400  samples.  The  choice  of  sequence 
length  is  related  to  the  specified  autocorrelation  function  and  spectral 
density  and  to  the  desired  resolution  in  time  and  frequency,  as 
explained  below. 

AUTOCORRELATION,  SPECTRAL  DENSITY  AND  RESOLUTION 

A  Gaussian  correlation  function  was  chosen  primarily  because  radar 
clutter  autocorrelation  functions  are  often  well  approximated  by  this 
shape  f 5 J .  Specifically  the  autocorrelation  function  was: 
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(21) 


R(t) 


.222 

-2ir  a  t 


with  corresponding  power  spectral  density  S(o>)  whose  shape  Is  also 
Gaussian: 


S(u>) 


ot/2rT 


/Q  2  2 
-u)/8n  a 


(22) 


Note  that  R(x)  is  normalized  such  that  R(0)  -  1  . 

A  At  of  1  was  chosen  for  the  correlation  function  samples  so 
that  Tn  *  nAx  =*  n  .  The  half  width  of  R(x )  is  defined  as  that  value 
of  x  for  which  R(x)/R(0)  -  .5  .  Clearly,  the  choice  of  x  depends 
upon  a  .  In  both  examples  the  value  of  a  was  chosen  so  that  the 
halfwidth  occurred  at  t  -  3  .  (This  is  in  rough  agreement  with  the 

number  of  pulses  over  which  clutter  echoes  are  typically  correlated) . 

2  2 

The  resulting  value  of  a  is  then  determined  by  a  »  in2/18n  .  The 

resolution  in  x  is  such  that  four  samples  (at  x  -  0,1, 2, 3)  fall 
within  the  halfwldth. 

The  half width  of  S(u>)  also  depends  on  a  ,  and  for  the  value 
specified  above  it  occurs  at  u>  »  (2in2)/3  .  Now  the  frequency  bin 
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width  Aw  is  equal  to  2x/NAt  ■  2n/N  ,  since  At  •  1  .  To  obtain  about 
the  same  resolution  in  w  as  in  t  one  merely  sets: 

m1/2  ,  T 1 / 2 
Aw  At 

where  Wj/2  »  Tl/2  are  tlie  halfw*dths  of  S(w)  and  R(t)  and  Aw  , 

At  the  respective  bin  widths.  Substituting  At  *  l  ,  Aw  *>  2it/N  and 
using  the  values  of  Wj/2  an<*  ri/2  mentioned  in  the  preceding 
paragraph,  a  sequence  length  N  a  41  samples  is  obtained.  Since  the 
Fast  Fourier  Transform  (FFT)  algorithm  used  to  compute  the  DFT  operates 
most  efficiently  when  N  is  a  power  of  2  ,  the  sequence  length  for 
this  analysis  was  chosen  to  be  the  next  higher  power  of  2  greater  than 
41  or  25  «  64  . 

Approximately  99.7%  of  the  total  power  is  contained  within  radian 
frequencies  less  than  or  equal  to  3  standard  deviations  of  the  Gaussian 
In  (22).  This  equates  to  an  actual  frequency  f  *  w/2x  of  .187  Hz. 
Twice  this  value  or  .375  Hz  would  be  the  corresponding  Nyquist 
frequency.  But  the  actual  sampling  frequency  is  1  Hz  so  the  spectrum  is 
essentially  free  from  aliasing. 

EXAMPLE  1,  SPECIFIED  AUTOCORRELATION  FUNCTION 

This  example  illustrates  the  generation  of  a  mean  1,  variance  1 
process  with  a  Gaussian  shaped  correlation  coefficient  as  expressed  in 


(21).  The  autocorrelation  function  differs  from  this  coefficient  only 
by  the  addition  of  the  square  of  the  mean  value  of  the  process.  The 
correlation  coefficient  is  shown  in  figure  la.  For  this  example  the 
uncorrelated  sequences  (XjJ  and  consequently  {U^}  were  generated 
according  to  a  uniform  probability  amplitude  distribution.  Samples  of 
the  spectral  density  given  in  (22)  and  shown  in  figure  lb,  were  used  to 
weight  the  .  A  DFT  was  performed  on  the  weighted  samples  to  obtain 
a  correlated,  mean  zero  time  sequence.  Finally,  each  sample  was  shifted 
by  one  unit  to  obtain  a  sequence  with  mean  1. 

The  sample  autocorrelation  function  (solid  curve)  is  compared  with 
the  desired  autocorrelation  function  (dashed  curve)  in  figure  2.  This 
sample  autocorrelation  was  computed  in  the  following  way.  For  each 
sequence  k  of  the  ensemble  an  estimate  of  the  correlation  function  at 
lag  l  was  made: 

*  k(i)  "ThT  Vi-K  (23> 

This  estimate  is  unbiased.  The  estimates  for  each  sequence  were  then 
averaged  over  the  ensemble  at  each  i: 

100 

‘  Too  k2l  Rk(i)  (24) 
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(a)  Correlation  coefficient  vs.  lag  r 


FIG.  1:  SPECIFIED  CORRELATION  FUNCTIONS  AND 
SPECTRAL  DENSITY  FOR  EXAMPLE  1 
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The  maximum  absolute  difference  between  these  curves  is  about  .05  or  5%. 

The  probability  density  function  is  shown  In  figure  3.  The 
histogram  is  a  frequency  plot  of  the  random  samples  and  the  solid  curve 
is  a  theoretical  normal  density  function  with  mean  and  standard 
deviation  equal  to  1.  This  theoretical  curve  is  normalized  to  the  area 
of  the  histogram.  The  left  hand  ordinate  scale  is  calibrated  in  numbers 
of  samples.  The  right  scale  gives  p(y)Ay  ,  the  theoretical  probability 
of  finding  a  sample  value  between  y  and  y-fAy  .  Here,  Ay  is  the  bin 
width  of  the  histogram.  Although  the  frequency  samples  were  uniformly 
distributed,  the  samples  in  figure  3  were  obtained  by  linear 
combinations  of  them  (the  DFT)  and  hence  by  the  Central  Limit  Theorem 
would  be  expected  to  tend  to  normality.  For  sequences  of  length  64,  the 
figure  indicates  that  agreement  with  the  normal  distribution  is 
excellent.  However,  as  clutter  samples  these  would  fail  in  some 
respects.  Although  they  are  appropriately  correlated  as  shown  in  figure 
2,  the  probability  density  is  not  characteristic  of  clutter  echoes.  Not 
only  do  16  percent  of  the  samples  have  negative  values,  but  the  long 
positive  tall  often  seen  in  clutter  distributions  is  not  present. 

EXAMPLE  2,  LOG-NORMAL  SEQUENCE  WITH  SPECIFIED  AUTOCORRELATION 

This  example  illustrates  the  generation  of  a  mean  1,  variance  1 
log-normal  process  with  the  same  Gaussian  shaped  correlation 
coefficient  Py(t)  as  used  for  example  1.  In  this  case  however,  the 
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spectral  density  of  the  underlying  normal  process  must  be  used  to  weight 

the  frequency  samples.  It  is  obtained  from  the  autocorrelation  function 

of  that  process,  R_(t)  which  is  related  to  p„(t)  by  equation  (4e) . 
z  y 

For  a  mean  1,  variance  1  log-normal  process  and  the  specific  form  of 
Py  used  here  (4e)  becomes 

T  .  2  2  2  "] 

R  (t)  -  in  [l  +  e  "2n  a  T  J  (25) 

z 


R  (t)  (solid  curve)  is  compared  with  p  (t)  (dashed  curve)  in  figure 
z  y 

4a.  Note  that  the  two  become  equal  for  r  greater  than  about  4.  This 
can  be  verified  quantitatively  by  expanding  ln(l+x)  for  small  x  . 


The  approximate  spectral  density  corresponding  to  (25)  wa6  obtained 
by  a  DFT  of  that  equation.  It  is  shown  in  figure  4b.  The  independent 
frequency  samples  and  {U^}  were  generated  according  to  a 

mean  zero  normal  distribution  and  weighted  with  appropriate  values  of 
this  spectral  density.  A  DFT  of  the  weighted  samples  resulted  in  a  mean 
zero  normally  distributed  time  sequence  whose  autocorrelation  function 
was  the  Rz(t)  of  equation  (25).  The  transformation 


U+1) 

y  *  e 


(26) 
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FIG.  4:  SPECIFIED  CORRELATION  FUNCTIONS  AND 
SPECTRAL  DENSITY  FOR  EXAMPLE  2 
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applied  Co  this  sequence  resulted  In  the  desired  log-normal  sequence 
with  mean  1,  variance  1  and  autocorrelation  function  Ry(t)  -  Py(t)+1  . 


The  sample  autocorrelation  function  (solid  curves)  is  compared  with 
the  specified  autocorrelation  (dashed  curve)  in  figure  5.  Two  ensembles 
of  100  sequences  each,  and  labeled  A  and  B  In  the  figure,  are  shown,  to 
illustrate  the  statistical  variability.  In  this  case,  the  maximum 
deviation  from  the  specified  function  is  about  13  percent.  The 
probability  distributions  of  each  of  these  ensembles  is  presented  In 
figures  6  and  7.  As  before  the  histograms  are  frequency  plots  of  the 
random  samples.  "Hie  solid  curves  are  theoretical  log-normal  density 
functions  with  mean  and  variance  both  equal  to  1.  The  scales  on  the 
left  and  right  hand  ordinates  are  exactly  as  described  In  figure  3.  The 
agreement  between  the  histograms  and  theoretical  curves  is  excellent  in 
both  cases. 


i:  SAMPLE  AUTOCORRELATION  FUNCTIONS  COMPARED  WITH  SPECIFIED 
AUTO  CORRELATION  FUNCTION  FOR  EXAMPLE  2 


(SAMPLE  A) 


.936 


probability  density  for  mean  1 
log-normal  process 

(SAMPLE  B) 
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